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ABSTRACT 



Aims. We examine the influence of nebular continuous and line emission in high redshift star forming galaxies on determinations of 
their age, formation redshift and other properties from SED fits. 

Methods. We include nebular emission consistently with the stellar emission in our SED fitting tool and analyse differentially a sample 
of 10 z ~ 6 galaxies in the GOODS-S field studied earlier by Eyles et al. (2007). 

Results. We find that the apparent Balmer/4000 A breaks observed in a number of z ~ 6 galaxies detected at ^ 3.6 jim with 
IRAC/Spitzer can be mimicked by the presence of strong restframe optical emission lines, implying in particular younger ages than 
previously thought. Applying these models to the small sample of z ~ 6 galaxies, we find that this effect may lead to a typical 
downward revision of their stellar ages by a factor ~ 3. In consequence their average formation redshift may drastically be reduced, 
and these objects may not have contributed to cosmic reionisation at z > 6. Extinction and stellar mass estimates may also be somewhat 
modified, but to a lesser extent. 

Conclusions. Careful SED fits including nebular emission and treating properly uncertainties and degeneracies are necessary for more 
accurate determinations of the physical parameters of high-z galaxies. 

Key words. Galaxies: starburst - Galaxies: ISM - Galaxies: high-redshift - Ultraviolet: galaxies 



1. Introduction 

Although star-forming high redshift galaxies have been identi- 
fied out to redshift z ^ 6, quite little is known about the phys- 
ical properties of these galaxies, since this requires in particu- 
lar sensitive observations in the near-IR range and beyond. With 
the advent of Spitzer it has recently become possible to detect 
galaxies at z ~ 6-7 at wavelengths ^3.6 //m, albeit still in quite 
small numb ers ( Egami et alJl2Q05b lEyles et ai]l2005l : lYan et alJ 
l2005l I2006I: lEvles et alJl2007h . One of main surprises of these 
studies has been the finding of a substantial Balmer/4000 A 
spectral break in a large fraction of these objects, taken as an 
indication for the presence of an "old" (several hundred Myr) 
underlying stellar population, which would indicate very high 
formation redsh ifts for these galaxies. For example, at z ~ 6, 
lYan et all d2006l) estimate representative ages of 40-500 Myr for 
their sample of ~ 50 IRAC detected objects; Hvles et alJ d2007l) 
analyse 10 IRAC detected galaxies and derive ages of ~ 200-700 
Myr and corresponding 7 < Zf orm < 18. However, it is possible 
that the Balmer break is affected or mimicked by emission lines 
in high-z galaxies with intense star formation, since their con- 
tribution to photometri c filters increases with (1 + z) (see e.g. 
Schaerer & Pello 2005). Determining properties such as maxi- 
mum ages, past star formation histories, and formation redshifts 
is important for a variety of topics including our knowledge of 
galaxy formation, our understanding of cosmic reionisation and 
its sources, and others. It is therefore essential to examine thor- 
oughly the uncertainties in the underlying SED analysis. 

While possible uncertainties in stellar mass estimates of 
high- z galaxies have the focus of several recent discussion s (see 
e.g. lMaraston et al.ll200o1; lBrazuaJl2007t lElsner et alj|2008l) . lit- 
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tie attention has been paid to the impact of nebular emission 
(both lines and continuous emission) on SED analysis, since 
these processes are usually not included in evolutionary syn- 
thesis models. The effect of nebular emission on broad-band 
photometry of galaxies has been studied in several papers (e.g . 
lAnders & Fritze-v. Alvenslebenl 120031; IZackrisson et alJ [2008). 
and it leads to i mprovements for photometric redshifts, a s 
shown recently by lllbertl d2009l) and iKotulla & Fritzel d2009l) . 
Nevertheless, such models have so far not b een applied to anal - 
yse the properties of high-z galaxies (but see lRaiter et al.l l2009). 
In this Letter we examine in a differential manner the effect of 
nebular emission on z ~ 6 galaxies, and we show that it can have 
a significant impact on determinations of their age and hence 
formation redshift. 

In Sect. [2] we summarise the galaxy sample and the SED 
fitting method. Our results and their implications are presented 
in Sect. [3] Our main conclusions are discussed and summarised 
in Sect. [4] We assume a flat ACDM cosmology with Hq = 70 
km s _1 Mpc" 1 , Q. M = 0.3, and Q vac = 0.7. 



2. Observational data and modeling tools 

2.1. Selection ofz~6 galaxies 

To examine the robustness of stellar ages and other physical pa- 
rameters of the most distant galaxies, we chose the sample of ten 
z ~ 6 star-for ming galaxies from the GOODS-South field mod- 
eled earlier by Eyle s et a"fl d2007l hereafter E07). Four of these 
objects have spectroscopic redshifts determined from their Lya 
emission; for the remaining objects we use the photometric red- 
shifts from the GOODS-MUSIC catalogue adopted by E07. We 
use the i'z'JK s , 3.6 //m, and 4.5 yum ACS, ISAAC, and IRAC 
photometry from E07. 
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Since longer wavelength data is also available, we have also 
used for comparison the IRAC photo metry at 5.8 and 8. fim 
from the GOODS-MUSIC catalogue of lGrazian etafl d2006l) for 
the subsample of 3 objects with spectroscopic redshifts also cov- 
ered by their data. 

2.2. SED fitting tool 

To analyse the SEDs we use a modified version of t he Hyperz 
photometr ic redshift code of Bolzonel la et all (120001) desc ribed 
e.g. in ISchaerer & Pellol d2005l) : ISchaerer et all (120071) and 
adapted to include nebular emission. Among the large diver- 
sity of spectral templates included in this v ersion, we here use 
the 20 03 GALAXEV synthesis models from lBruzual & Chariot! 
(120031) . covering different metallicities and a wide range of star 
formation (SF) histories (bursts, exponentially decreasing, or 
constant SF). For comparison with E07, we define the stellar 
age f* as the age since the onset of star formation, we adopt a 
Salpeter IMF from 0.1 to 100 M Q , and we properly treat the re- 
turned ISM mass from stars. 

To account for the effects of nebular emission from young, 
massive stars on the SED we include nebular emission (both 
lines and continua) in a simple manner. Continuum emission 
is added to the stellar SEP, as in our synthesis models (cf . 
ISchaerer & Vaccalll998t lEeitherer et al]fl999t ISchaereril2003h . 
The main emission lines of He, C, N, O, S, and other lines 
are inclu ded using the empirical relative line i ntensities com- 
piled by lAnders & Fritze-v. Alvenslebenl d2003l) from galaxies 
grouped in three metallicity intervals covering ~ 1/50 Z to Z Q Q 
In addition we include H recombination lines from the Balmer, 
Paschen, and Brackett s eries, as well as Lyq; the ir relative in- 
tensities are taken from Stor ey & Hummer! (1 19951) for a typical 
ISM (density,temperature) of (n e = 100 cm" 3 , 10 4 K).0 The ab- 
solute strength of both the continuous and line emission depend 
to first order only on the total number of Lyman continuum pho- 
tons, which can be computed from the template stellar SED. In 
this manner we include the main nebular emission features from 
the UV (Lya) to 2 yum (restframe), necessary to fits the SED of 
galaxies at z > 4 up to 10 //m (IRAC Channel 4). 

The free parameters of our SED fits are: the metallicity Z 
(of stars and gas), the SF history described by the timescale t 
(where the SF rate is SFR oc exp ^ T ), the age t, the extin ction 
Ay described here by the Calzetti law ( Calzetti et al.ll2000l) . and 
whether or not nebular emission is included. Here we consider 
three metallicities Z/Z Q =1, 1/5, 1/20, a wide range of t values 
as well as bursts and SFR=const, ages up to the Hubble time, 
and Ay — 0-2 mag. More details on th e SED fitting method will 
be presented in lde Barros et a i1 d2009l) . 



3. Results 

3. 1 . Comparison of SED fits with/without nebular emission 

Overall we obtain similar although not identical results to E07, 
when using the same assumptions. In other words, adopting solar 
metallicity and neglecting nebular emission we find that a signif- 
icant fraction of the objects (at least 5/10) are best fit by fairly 
large ages of the order of f* <: 500 Myr and generally very low 
extinction. The quality of our fits (expressed in reduced x 2 ) is 
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Fig. 1. Observed (blue points) and best fit SEDs (solid lines) 
of the z = 5.83 galaxy 23_ 6714 from E07. The errorbars of 
the observed wavelength indicate the surface of the normalised 
filter transmission curve. Upper limits in flux indicate l<x lim- 
its. Red crosses show the synthesised flux in the filters. The two 
SED fits shown are based on a standard Bruzual & Chariot so- 
lar metallicity model (magenta), and the same modeling includ- 
ing also nebular emission (black). While the age of the former 
(with^ 2 = 1.61 ) is £ 700 Myr, the latter gives U ~ 20 Myr 
Or 2 =0.14). 



comparable to E07, and similar values are also found for the stel- 
lar masses and current SFR. For example, the mean age, stellar 
mass, and extinction of the 10 objects we obtain with these as- 
sumptions are U ~ 500 Myr, = 1.2xl0 10 M o , and Ay' = 0.11 
(i.e. I^v ~ 0.05). 

However, including nebular emission changes quite drasti- 
cally these results. In this case only 1 of 10 objects has a best fit 
age ti, ^ 500 Myr, and the average age of the sample is lowered 
by a factor ~ 4 to U ~ 120 Myr. A slightly lower mean stellar 
mass M* = 7.8 x 10 9 M , and a higher extinction Ay = 0.34 
Eb-v ~ 0.08) is also obtained. Compared to the above set of 
models the^- 2 values is found to be lower for half of the objects. 
A more rigorous comparison, including also a careful discussion 
of the uncertainties, is deferred to a later paper. 

Examples of SED fits for two objects with known Lya emis- 
sion and known spectroscopic redshift are show in Figs. 1 to 3. 
Figure Q] shows the best fit solution using pure stellar SEDs at 
solar metallicity compared to the best fit including nebular emis- 
sion (found to be for Z = 1/5 Z ). While the age of the latter 
is t+ « 20 Myr, the former has an age close to the maximum 
allowed for the redshift of this object (z - 5.8). With nebular 
emission the apparent Balmer break is mimicked by the pres- 
ence of strong restframe optical emission lines boosting the flux 
both in the 3.6 and the 4.5 fjm filters^ For this object (23-6714) 



1 See lKotulla etai] J2009I) for a comparison of the resulting emission 
line spectrum with nearby galaxies. 

2 The emissivities of these lines are known to depend little on n e and 
T. 



3 The strongest lines in the 3.6 fim filter are [O m] /1/14959,5007 and 
H/3, the strongest at ~ 4.5 yimi are Ha, [NH]/l/l6548,6584, and [S ill] 
/L19069,9532. 
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Fig. 2. Same as Fig. 1 for 23.6714, but including IRAC 5.8 and 
8.0 yum data (l<x limits), and using the GOODS-MUSIC photom- 
etry. The magenta line shows the solar metallicity "standard" 
SED fit with an old age (f* ^ 700 Myr) and no extinction ob- 
tained in Fig. 1. The black line shows our best fit including neb- 
ular emission with a very young age (f* ~ 6 Myr) and Ay ~ 0.7, 
the green line the <: 700 Myr best fit with standard (purely stel- 
lar) templates. 



Fig.3.Same as Fig. 1 for object 31J2185. The two small blue cir- 
cles illustrated the difference in the photometry of E07 compared 
to that of the GOODS-MUSIC catalogue shown with large sym- 
bols. The red (magenta) line shows the best fit to the photometry 
of E07 assuming solar metallicity models with (without) nebular 
emission. The black (green) lines show best fits to the GOODS- 
MUSIC photometry with (without) nebular emission. For both 
datasets the age obtained with nebular emission is f* ~ 90-180 
Myr, compared to ~ 130-700 Myr with standard templates. 



a non-zero extinction, of the order of Ay ~ 0.2-0.8 depending 
on metallicity, is required. 

Since nebular lines at this redshift are expected to affect less 
the photometry at longer wavelenghts (^ 5 /mi), it is interesting 
to examine SED fits including also the available 5.8 and 8.0 /urn 
data (see Fig.|2|. For consistency we therefore use the GOODS- 
MUSIC photometry in all bands. A good agreement with the 
observations is found, and again significantly lower ages are ob- 
tained with nebular emission. For 23-6714 we find the best fit 
with Z = 1/5 Z Q , a very young age of f* ~ 6 Myr and Ay ~ 0.7, 
compared to ages of ^ 700 Myr with standard (purely stellar) 
templates. 

Figure[3]shows an object (31-2185 from E07) with an appar- 
ent Balmer break of similar strength. An old age f* ~ 700 Myr 
is again obtained here for standard stellar SED fits, whereas the 
inclusion of nebular emission lowers the age to 90-180 Myr, de- 
pending if we include or not the photometry in Channels 3 and 4. 
For this object and for the GOODS-MUSIC photometry it turns 
out that nebular emission is not crucial, since relatively short 
star-formation time scales and hence solutions with few ionising 
stars are preferred. 

The above results clearly illustrate that it is important to in- 
clude the effects of nebular emission (lines in particular) to de- 
termine the properties of these high-z galaxies, most importantly 
their stellar ages. 

3.1 .1 . Comparison with earlier work 



are crucial to measure the Balmer break or its absence? There are 
three main reasons for it. First, E07 compute the contribution 
of Ha to the 4.5 yum filter from an estimate of the unobscured 
SFR(UV), which assumes a long (S 100 Myr) SF timescale. 
However, for shorter timescales the Ha lumino sity can be larger 
by up to a factor ~ 4, as illustrated e.g. by (IVerhamme et alj 
2008, see their Fig. 15). Second, as noted by E07 the contribution 
of Ha is further enhanced when extinction is taken into account. 
Finally, [O in] /1/14959,5007 plus H/3 have strong enough intrin- 
sic intensities (typically 1.8-2.2 times Ha) to boost as well the 
flux in the 3.6 fim filter, thereby mimicking a "Balmer break" be- 
tween IRAC Channels 1-2 and shorter wavelength filters. When 
the effects of nebular emission lines are taken into account in a 
consistent manner with the ionising flux predicted by the SF his- 
tory and when extinction is allowed for, we find that - at z ~ 6 
considered here - their contribution to IRAC channels 1 and 2 
can be significant This largely explains our lower age estimates 
compared to the earlier work of E07. 

3.1 .2. How to distinguish "old" and young populations? 

Can the two type of SED fits, explaining the Balmer break with 
"old" stars or by the presence of emission lines, be distinguished 
observationally? 

For the SED fit of 23-6714 shown in Fig.[TJfor example, the 
current SFR ~ 800 M Q yr _1 . From this the expected Ha flux is 
of the order 3. x 10~ 16 erg s _1 cirT 2 , clearly beyond the reach 



Why does this conclusion disagree with E07 who also examined 
the possibility of "line contamination" in the IRAC filters, which 



4 See also the predictions of Zackr isson et al.1 d2008l) who also find 
significant line contamination even for constant SF over 50 Myr. 
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of current spectroscopic facilities at~ 4.5 fim. This object being 
the brightest of the sample it appears that direct detections of the 
restframe optical emission lines need to await future facilities, 
such as the JWST. 

In principle deep photometry at longer wavelengths, where 
the expected strength of emission lines decreases, could help to 
disentangle the two solutions. However, the current 5.8 and 8.0 
/im data from Spitzer is not deep enough to rule out the "old" 
SEDs, as shown for the two brightest objects 23-6714 and 31- 
2185. Although the spectrum of the ~ 700 Myr old popula- 
tion displays an excess with respect to the 5.8 and/or 8.0 fim 
fluxes, its significance is 5= 2<x at best. Future, deeper observa- 
tions should be able to set firmer limits on this issue. 

In any case we note that neglecting nebular emission for 
high-z galaxies with ongoing massive star formation — as testi- 
fied e.g. by their Lya emission often used to confirm their red- 
shift spectroscopically — is physically inconsistent and its ef- 
fect on the determination of their physical parameters should be 
taken into account. 

3.2. Implications 

As shown above, the inclusion of nebular emission (lines and 
continuum) can alter the physical properties of galaxies deter- 
mined from their SED. For example, considering just the average 
properties of our best fit solutions for the 10 objects fitted here 
we obtain U ~ 400 Myr, = 1 .1 x 10 10 M G , and = 0. 16 for 
standard Bruzual & Chariot spectral templates when metallicity 
is also varied. Including nebular emission we obtain f+ » 120 
Myr, M7 = 7.9 x 10 9 M , and a higher extinction AjJ = 0.34. 
Although indicative of the trend obtained with nebular emission, 
the reader should be aware that large deviations are obtained and 
uncertainties should be treated properly (cf. Sect. |4j. 

Clearly the most striking result is that the average age of the 
z ~ 6 galaxies may be decr eased typically by a factor ~ 3 com - 
pared to earlier studies (cf. lYan et aT]|2006t Hvles et al.ll2007l) . 
Translated to the average redshift of formation Zf orm this would 
imply a shift from Zf orm ~ 9. to 6.6, assuming z = 5.873 for the 
average redshift of our sample (cf. E07). For only one out of 10 
objects do we find best fit ages of > 200 Myr, i.e. a formation 
redshift beyond 7. If true, this means in particular that the con- 
tribution of the galaxies currently observed at z ~ 6 to cosmic 
reionisation must have been negligible at z ^ 7, in contrast to 
the finding of E07. 

The average extinction of Ay ~ 0.34 mag (£b-v ~ 0.08) 
we find, corresponding to a UV attenuation by a factor ~ 2.2 
is also worth noticing. With standard spectral templates, ne- 
glecting nebular emission, the average is Ay = 0.16. This re- 
sult indicates that dust attenuation may be larger than previously 
thou ght for Lyman Break galaxies at this redshift. For compari- 
son, iBouwens et alj d2008l) assume a UV attenuation factor of ~ 
1.5 at z ~ 6. If true and representative for the population of z ~ 6 
galaxies, it may imply an upward revision of the SFR density 
by ~ 50 %. Finally stellar masses and hence the estimated stel- 
lar mass density of i-drop galaxies may also be slightly affected. 
Also already mentioned, the average stellar mass of our sample 
Mi, is ~ 30 % lower when we include nebular emission in the 
SED fits and leave the other assumptions unchanged. 

4. Discussion and conclusions 

The main objective of this work has been to examine the effect of 
nebular emission (lines and continua) on the derivation of phys- 



ical parameters of z ~ 6 galaxies through broad band SED fits. 
Although our work clearly indicates that this effect can signif- 
icantly alter derived properties such as galaxy ages and extinc- 
tion, we do not (yet) aim at determining absolute values. 

Indeed both uncertainties/difficulties in the photometry as 
well as uncertainties and degeneracies in the models remain and 
their effect needs to be quantified properly. For example, the dif- 
ferences obtained using photometry from two sources (revisions) 
has been illustrated above. Stellar ages and other quantities may 
also vary if different, more complex SF histories, different ex- 
tinction laws (possibly also extinction differences between stars 
and the gas), different IMFs are allowed in the models. 

Last but not least, our and other SED fitting tools are also 
limited by their dependence on stellar tracks, which are e.g. 
sensitive to input p hysics such as the treatment of rotation 
dVazquez et alj|2007b. to the treatment of advanced phases such 
as TP-AGB stars dMaraston et al.l20 06, cf.) and others. Once the 
input physics of the SED fitting code and the photometric data is 
fixed, errors and degeneracies in the fit parameters need to quan- 
tified in detail and for each individual object. Finally appropriate 
quantitative methods should be used to determine the average 
properties and their uncertainties for samples of galaxies. Such 
detailed analysis, clearly beyond the sc ope of this Letter, will be 
presented in a forthcoming publication dde Barros et al . 2009). 

Independently of these remaining uncertainties, our work 
shows clearly that the presence of strong B aimer/4000 A breaks 
in z ~ 6 galaxies is not necessarily due to stars, but can 
also be mimicked by restframe optical nebular emission lines, 
whose contribution to broad band filters increases with (1 + z). 
Apparently "old" ages and high formation redshifts of galaxies 
with strong ongoing star formation may need to be revised sig- 
nificantly downward. Finally, the impact of nebular emission on 
SED analysis of other galaxies, also at lower redshift, needs also 
to be examined systematically. 
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